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Description. A classical problem in physics and geometry is the qualitative description of the spectrum of 

C*~) ■ the laplacian on functions on a domain with Dirichlet boundary conditions. The eigenfunctions determine the 

standing vibrational modes of a drum shaped as the given domain. It is expected that the eigenfunctions behave 

iT^ ■ as polynomials of degree related to the order of the harmonic. This is manifested for instance by the size of their 

nodal sets, i.e. the set of zeroes. In terms of the vibrational motion of a drum it consists of points that remain 

stationary and can be interpreted as the locus of destructive interference of the waves with the boundary. In 

the late 18th century Chladni performed experiments with planar drums of that revealed the shape and size of 

the nodal set vary with the order of the harmonic Similar questions can be asked for the laplacian acting on 

^~j | functions on a an arbitrary riemannian compact, closed manifold 

Statement and known results Let (M. n ,g),n > 2 be an n-dimensional compact closed manifold equipped 
with a smooth riemannian metric g. The Laplace-Beltrami operator acting on functions A g on M n is written 
"*^ in local coordinates, g l i = (g~ 1 ) i j,g = det(gij): 



^ 



> 



X 



g: ^-^IM ! < 



The spectrum of the laplacian is discrete spec(A) = {Aj}., 6 n C (0, oo) consists of eigenvalues with eigenfunctions 
ma, Ae spec(A) are the solutions of the equation 

Agu x + Xu x = 



The nodal set of an eigenfunction u\ is defined as 

l> ■ N(« A ) = {x G M n /u x (x) = 0} 

and it is proved that this set is a smooth submanifold outside a set of (n— l)-Hausdorff measure zero. Briinning 
£f~) ■ [1979] and [Y] showed that for a smooth surface the nodal length has a lower bound as: 



jr x (N(« A )) > CA 1/a 

and then Yau conjectured (|Y|. pr. no 74) that the Hausdorff measure of this set grows as 



H ■ jr n - 1 (N(u A )) ~ A 1/2 

c3 ■ 

which is evidently the case for the spherical harmonics on the sphere. [DFlj established this when (M. n ,g) is 
analytic. Furthermore they showed [DF2 that for smooth surfaces one has that 

J?\N(ux)) <CX 3/4 

which was obtained by ! D by different methods. Note that [N] proved that 

^(NOa)) <C"AlogA 

Hardt and Simon in their excellent work |HS] showed for C 1 ' 1 metrics that 

JM based on Donnely-Fefferman work obtained bounds for the size of tubular neighbouhoods of nodal sets for 
real analytic metrics. Recently |SZ) . [CM] . [M] came up with new lower estimates of the nodal volume. The 
first authors provide an interesting formula for the size of level sets comprising an integral over the regular level 
sets of the eigenfunction. 



Short description of the method. Our method is based on the construction of a cell decomposition of a 
riemannian manifold with the boundaries of cells consisting of pieces of geodesic spheres. This resembles the 
picture suggested by the Huyghens principle governing wave propagation. The cells are selected so that we 
are allowed to follow an inductive argument restricting the eigenfunction on the faces and get precise growth 
estimates for it. Recall that the area of geodesic spheres determines the absorption rate of a wave in the course 
of propagation in a medium. The faces are constructed through the introduction of a cluster of points: these 
are thought as a collection of sources of spherical waves, they define the geodesic spheres that provide the pieces 
forming the cells. The domains that are formed are called geodesic pixels while the faces are called undulating 
fronts and the geodesic sphere pieces called elementary wave fronts. The pixel size is arranged essentially by 
the mean curvature of the elementary wave front. The conjecture is then proved by an inductive argument 
combined in the lower bound with the isoperimetric inequality and eigenvalue estimates. The upper bound is 
obtained by an elaboration of the Dong formula. 

Notation We introduce some notation. First we introduce the localized energy of a geodesic pixel: 

E«(R,B;C)= f C 2 (I VR | 2 + I B | 2 ) 

and 

E<°>(R;0= / C 2 |R| 3 
where we introduce the Bach tensor 

Bl3k = ^—^ CUTl ( Ric )^k - 2 ( n -l)(n-2) ^ 9ijRk ~ ftfcR ^ 

The faces of the pixel are pieces of geodesic spheres called elementary wave fronts and are denoted by F with 
mean curvature h and localized tension T (/i;$) given for a smooth test function #, supp$ C F: 



2.2 



T 1 (h;&) = / $ 2 \^h\ 2 , T°(M)= / $ 2 h 



Then we introduce two numbers 



and 



C6C«(P)V E°(R;C) 
t(F) = sup ! 



*EC ~(F)VT°(M), 

The first theorem comprises Harnack estimates on such sets 

Theorem 1. Let (M",j) be a riemannian manifold with scalar curvature functional. Then in apixelP, r(P) = 
{£; i\, . . . , ik) =11 with generic face F, i(F) — r\ we have that for positive constants explicitly calculated 

inf \h\ < en inf \h\ + e, 

F F 

inf \k\ 2 < ci 2 inf \k\ 2 + e 

F F 

where Cn,ci2 depend in an explicit way on F. Furthermore let u be an eigenfunction with eigenvalue A and e a 
regular value of u 

P £ = Pn{|u A | >e}, F £ = Fn{|u A | >e} 

Then for explicitly calculated positive constants C2j — C2j(P, A), c^j = Csj(F),j = 0, 1, 2, 3, 4 we have that 

sup \u\ < C 2 o£ 
p £ 

sup(|Vu|) < C21 inf |Vw| + c 2 4 
p. p - 

sup(|V 2 u|) < c 22 inf |V 2 u| + c 25 
p £ p e 

sup|u| < c 30 e 

F t 

sup(|yu|) < c 3 i inf \yu\ + c 34 

F t F ^ 

sup(^ 2 u|) < c 32 inf \^u\ + c 35 

F. F ^ 



These estimates provide the basis for the inductive argument and the application of Dong's formula: 
Theorem 2. Let u be as above then 

1 Decomposition of a manifold 

1.1 Quantitative Huygens principle 

We introduce the set of points, chosen to lie on the nodal set 

Co = {<2°}£i c M" 
and lying at distance 

This is possible since in every geodesic ball of radius 0{\~?) there is always a zero. Furthrmorc we define 
geodesic balls of radii r^, B™ bounded by geodesic spheres S™^" 1 . These balls are taken to overlap in 
(n + l)-ples : we set as p°;*i ■■■*<= the overlap region of the balls : 



■pO;ii...ik _ r~\ dii 



3 = 1 

for k = 2, . . . ,n + 1. These domains are bounded by spherical regions denoted as F '* 1 ""'*' , and have interior 
curvature data 

and face second fundamental form 

i/0;ii...i k -,e gO;ii...i k ;£ 

We call these sets geodesic pixels. 

We produce generations of such pixels after the introduction of new centers and arrive at the collection of 

pixels after generation j: 

pj;*i -ifc 

Its boundary consists of the elementary wave fronts F j;n '-^ fc;£ and curvature data: 

Ric j ' il - ik , fin-*, p-M-iK^ 

It is written in the form for 

Qpj\ii...i k _ | I j?j';ii...i k ,;t' 
j',k',i' 

The faces F- 7 '^■■■ z k'< i are called elementary wave fronts (EWF). Each pixel defines homothctic EWF spanned 
by the tubular neighbouhoods of the elementary wave fronts: 

K (j;M)q = Irs x p/;ii...i» f j re = ((i _ £ y } (i + e ) r ) 
We introduce the localized tension T(/i;i?) of an EWF given for a smooth test function -#,supp-# C F ;n ' 4fc : 

V(h;tf) = f $\tf j h\ 2 



where fi is the mean curvature of the (EWF) . We introduce two numbers 

'E (1) (R,B;C)' 



r 



(P) = sup 



C6cF(P)V E0(R;C) 

and 

t(F) = sup ' 



eec o ~(F)\T°(/i;0) 



Let rig, /j, be positive constants. We say that a pixel P with boundary consisting of EWF 

n+1 

OP = (J F e 

e=i 

satisfies an (%,/tt) condition if 

t(F e ) < m, r(P) < M 

Let (rij' ; i 1 ...i k ,;£,^j, e e) be positive constants. We assume that the geodesic pixels pO';*i>— ,»/) are selected so 
that their EWF F^ i*i ■■■**' sO satify (j/j/;^...*,, ,q) estimate while (/1j,q) are the parameters in the curvature 
estimate. This set of pixels is a subset of compact closure in the neighbourhood of the zero section in TM. 
In the sequel we will try to establish the equations defining this set, as consequences of the above integral 
inequalities. 

The structure of the metric in geodesic polar coordinates. Every geodesic pixel has center some C? 
obtained in the j-th generation of center selection. We consider geodesic coordinates from the neighboring 
pixels. Therefore let B"^ be a geodesic ball centered at the point C? and introduce polar coordinates through 

Gauss lemma. The metric is written then as: 

g = dr 2 + jf(r) 

where i{r) is a riemannian metric on the geodesic sphere S"~ = <9B™ t with second fundamental form and 

mean curvature respectively /c, fi. Accordingly we have the first and second variation equations for the metric 7, 
if we denote the radial derivative by 4- while the angular ones by ; j while in this section the symbol Jv denotes 
collectively the angular derivatives: 

(la) 

(lb) 
ar 

Then setting: 



we infer that 



df) 
dr 


= 2k, 






(1 n j j 

dr 


— — k 


urn 
K 3 


' R-OiOj 


i 


= \\og 


(det( 7 )) 


a = R i0 jok ij + k l m 






K = Kq 


= 1*1 






d:fi 

dr 


= h 






dn 








dr 


= a 





(2) 



Based on Newton's identities for symmetric polynomials we get the following inequality for a: 
(3 - C n , 3 )hn 3 -h 3 - |Rm|« < o < Cn ' 3 + 3 hK 2 - \h 3 + |Rm|« 



or 

-(C„,3 + 2n)n 3 - \Rm\n < a < C "-3 + 3ra K 3 + | Rm | K 

O 

The structure equations The Gauss equations that relate the curvature of 7, R to the ambient curvature 

-K'imjn ~r~ \fcij fcmn ^in^jm ) — ^imjnj \0<i) 

R + h 2 -k 2 = R-R 00 , (3c) 

The Codazzi equations 



constitute a Hodge system: 

.CUTl^/C Jijni tXraoij , -. 

dJv(k)i - Vih = R Qi 
where 

1.2 The second fundamental form and the mean curvature of the fronts 
1.2.1 Harnack on the slice 

The Gauss-Codazzi equation for the spherical front gives that 

\yh\ 2 = jM{k)ihi - Roihi, 

(5) 



Elaborating the preceding identities with Young inequality for sutiable p and obtain /i-growth inequality in the 
spherical front domain F = F <J;n "- lfc:£ ) and cut-off #: 

f Wfe) | 2 < 2f f *>|* + (p-iM^m) 1 ^ /■ ^ W | 2 

Jf P Jf pep- 1 Jf 

or choosing e = ( ^ p ~ ' j P vol(J^fc^) 2 p 7 and get that 

/" ^(^)| 2 < T) (- J (p- lf-^ol^fc,^) 2 ^ / tf 2 fr 2 P 

We recall Sobolev inequality from |SLlj for the case of EWF, F C W, r = ""l^ 

(j^r)" <cjT iyj7| + |fc||j7| 

Starting from 

(7 IDT 6 ") " < «? /" \uf- l \yu\ + \h\\U\ e 

and applying Holder's inequality 



J \u\ i - i ^u\ < (voi(F)) 1 -^ (J iu^-A" (J \yu\^ 



i—i 



(6) 



where q < 2. The inequalities on the slice 

/ ^l 2 < C / \&f{U)\ 2 + b«^(^)| 2 + \Rm\\U\ 2 
Jf Jf 

read for the localization second fundamental form, U = £|fc|: 

c 2 W < c f \c\ 2 \m c \ + c\yh\ 2 + (\yc\ 2 + irhc 2 ) 1*1 

JF 



The last term in the right hand side leads after the application of Young's inequality combined with Sobolev's 
inequality to 

/ C 2 W < C f |CI 2 |R*c| + c^| 2 + / {WC\ 2 + IRM)^ 1 

Jf Jf Jf 

We consider now the slice regions Fj determined by the tension energy through the sequence of constants {t]j}, 
in which t(Fj) — rjj. The Harnack inequality is proved through Moser iteration on the domain with smooth 



boundary W C F obtained by smoothing out the boundary of Fj . Therefore we exhaust the domain through 
the harmonic approximation of the face defining function F, Fq : 

Wjfa) = {x e F : (9 - 9^ +l )ri < \%{x)\ < (1 - 6 + 6 j )r]} 

and then Harnack inequality takes the form: 

sup \h\ < D(r],r]e) inf \h\ 

where the quantity D(ji,r]e) > is calculated in the appendix. 

The growth of the tension integral The preceding estimates necessitate the derivation of radial growth 
estimates for the tension integral 

T = T 1 (h;d) 

We start differentiating and proceed with the application of the structural equations. We obtain the differentail 
inequality 



dr) 
dr 



We obtain that 



d ^<2(fd 2 \h\$fh\ 2 +f$fh\ 2 d 

<2|fc|^fc| + ^R 00 



+ / r\h i \ 

F 



dhi 



dr 



dhi 

dr 



Therefore we have that: 



dT 

dr 



< (sup \h\\ t + (sup m) Qf \^h\^ U ly^i 2 " 



The last term is majorised by 



X^' 2 ) (X^' 2 ) ^ c,t +X^ 2|r * c|+|rw| ^ c|2 



Furthermore we have that: 

We select cutoffs satisfying for some rj > 

We arrive at the inequality 



SUp |fc| < CiT 2 + C2 
F 



d\tfd\ 



dr 



<ri&#\ 



dT 

dr 



<ciT 3 / 2 + c 2 , Ci = Ci(F,e), t = l,2 



We conclude through the use of Young's inequality: 



and hence for r < JS,: 



dT < 9 s 

— < c (T 2 + 1) 
dr v ' 



|T(r(l + e)) - T(r(l - e))| < tan(CV) < Cr 



1.2.2 Radial Harnack estimates 

In this section we derive the radial variation of the curvature quantities in the radial interval J rjE = ((1 — e)r, (1 + 
e)r) relative to a given value of this quatities at r. We write (|lbj) in contracted form: 



dh 
dr 



= -k 2 - R, 



do 



(7) 



Fisrt we have that 



— < -h 2 - R 00 & h 2 < - n - 1 )— - (n - l)Roo 

dr n — 1 dr 



We derive the differential inequality for k 

d\k\ 



dr 



< |fcr + |Rm|<(|fc|^ + l)(|Rm| + l) 



that is written after integration and elementary trigonometry as: 

||fc(r(l +£ ))|-|fc( r (l_ £ ))|| 



< tan(Cr) 



for 



For r < 97r ■ 
rul ' - 20C- 



l + |fc|(r(l + e))|*|(r(l-e)) 
max(C(e) / |Rm| ] ,max(|fc(r(l + e))\ ■ \k(r(l - e)|) < C 



K 



\\k(r(l + e))\ - \k(r(l - e))\\ <Cr 



The estimates for |Vfc|, |V 2 fc| follow from the elementary differential inequalities: 



d|Vfc| 
dr 



< \k\\Vk\ + (|VRm| + |Rm||fc| 



rf|V 2 fc| 
dr 
This lead to the estimates: 



Similarly 



< 



C (|fc||V 2 fc| + |Vfc| 2 + |Rm||Vfc| + |VRm||fc| + |V 2 Rm|) 



|Vfc| < 



|V 2 fc| < 



Cr 2 | / |Rm| 



K 



IVRml 



K 



Cr 2 



IRtoI 



K 



IVRml 



K 



Let x be a cutoff supported in supp(x) C I r .e such that 

e|x'|+e 2 |x"l<C 

as well as 

Me (r)= / X h(r',Odr' 



Lemma 3. The following estimates hold for r < —h= and are relative to a fixed value of h(r): 

\^\<(2sr)^ 2 (\h\'/ 2 + (2er\C \y/ 2 

[VMel < 2er([y/i| + 2(l + e)rci(|Rm|,|VRm|)) 

I^Vl < (2er) (\^h\ + 2r(l + e)c 2 (|Rm|, |VRm|, |V 2 Rm|)) 

The estimate of /j, £ We start applying and Cauchy-Schwarz in the end 

iMsl < (2er) 1 / 2 I f h 2 \ <{{n-l)er) l l 2 {\h{r{l-e))-h{r{l+e))\ 1 l 2 + {2eC r) 1 / 2 ^ < (n-l)*re(|fc|+C ) 

Estimate of [Xa* £ |, \y^/i e \ We commence with the integration by parts for U : I r . e — >• R: 

d\U\ 



f 




f dr 


/ u 


= 


/ -rC/ 


Jzr,. 




h re dr 



and apply this to 



d\yh\ 



dr 



<) '%w\<*m-l 



<Ci(|fc|p/fc| + ^Roo|) 



dr 



dr 



<C 2 {\k\^k\ + ^k\ 2 + \^R m \) 



and obtain the desired estimates. 



The radial-slice estimates. We start recalling the differentiation identity: 



<(/„W? + M7 



dr 



dr 



'F / JF 

We introduce a cutoff function \ supported in the shell suppx C K: 

K = I r . e x F, 
I r , e = ((1 - e)r, (1 + e)r) 

Therefore applying Sobolev inequality on the slice, Holder and standard elliptic estimates we derive the differ- 
ential inequalities for the quantities 



—± <aU. +bj, 
dr 



dH 



-r- < aJU~H j + c, 
dr v j j 



where 



U, 



/"-I ! 



jfc| S 



Hi 



iy^i : 



a = D{rj) 



\k\ 2 + h 2 



-,1/2 



M 



6i = / \VRm\ z + U / |Rm| 
'f Jf 



Ui + Uo I |Rm| 

F 



62 = ^1 

The diffrential equations for Uj , -ffj are of the form 



|Rm|, 



IVRml 2 U< 



H) 



|V 2 Rm| 2 



(8a) 

(8b) 
(8c) 
(8d) 



dU TT , 

— <aU + b 
dr 



we obtain for \g\ < e: 



U(r(l + g)) < e r ti a ( 1 + r "> dr ' ( U(r) +r 6(1 + r')e- r ' SS' ^+e'W 



Lemma 4. The following hold true in I r£ for i,j = 1,2 and constants dj = Cij(r; 771, ... , r/j-; e.,) > 

17 if ±(r(l ± e)) < c^.^r), H jt± (r(l ± e)) < C2jH j>± (r) 

1.3 Selection of the pixels through the front estimates 

1.3.1 The basic ansatz 

We introduce the notation for a smooth domain W in a riemannian manifold equipped with riemannian volume 
dv: 

V j (U:W) = J \V j U\ 2 dv 

We will derive the equation satisfied by u\ near (EWF) i.e. in a spherical shell of the form 

K = I r%E x F = ((1 - e)r, (1 + e)r) x F 

for some F = F^' 11 ""'*'^ with coordinates (r,6) and volume drda — y/jdrd6, while hi stands for the mean 
curvature of the shell. 

We select the front so that its mean curvature is controlled by the eigcnfunction growth near it. The 
eigenfunction equation is written near (EWF) in the form: 



d 2 u x du\ . A 

—tt- + h— h A 7 m a 

dr A dr 



Xu\ 



We make the ansatz for the parameter /3 to be determined and the smooth functions a, < 

u x (r,6) = A (r, 6) sin (M(r,e)) , 



and arrange the (EWF) so that near the front A is smooth and positive which implies that <p inherits the unique 
continuation property from u\. The equation then splits as: 



d 2 A h dA_ 2 

dr 2 dr \ V dr 



^>| 2 )A + (X-A + XA) 



sin(/3i0)- 



+Pi 



2^ + hA )^+A^ + (2^A-^ + A^) 
dr I dr dr z 



cos(/3i0) = 



This suggests after the selection of suitable positive function a depending on the angular variables: 



and consequently: 



A(V\r, 0) = e A , A = -~ Me + p 2 a(6) 



^(^) 2 + ^M 2 = iK + A-XA-|yA| 2 



d 2 ^ 
dr 2 



-yA-yv> = o 



where 



K = h 2 + 2(\k\ 2 + \R 00 \) 



(9) 



(10a) 
(10b) 



Then (jlOap provides the pointwise estimate: 

- /3 2 aX« ~ (1 + e)f3la\tfa\ 2 < (|K| + A) a - (j^x + \W^A a (11) 

We employ Hardy's inequalities for the harmonic approximation of p, in the terms of the last parentheses: 



J \£n + \^n » < sup \n\ J 

provided that 



Am 


a + -(sup|^|) 2 / 

^ F Jf 

SUp |/i| < £Q 
F 





J — a < e (c 3 (H) + i C2 (iJ) sup | M | J J tfo 



Set also 



The second equation gives that 



f d 2 \^a\ 2 <4 / (|K|+A)tf 2 

F JF 



1o = e 2 + f3 2 e 2 2Ai = LV>/| 2 + %^a| 2 



i^A^ + ^l 2 



(12) 



(13) 



We introduce arbitrary test functions in the angular variables: ■& supported in F. Then we multiply the equation 
by t} 2 and integrate using the identities: 



u d -^ = li([uA-llhu* 



F dr 2 dr \Jp ) 2 j F 



dU\ z d 



dh 

dr 



— = -\k\ 2 -R, 



0(1 



(14a) 
(14b) 
(14c) 



Hence we obtain that 



dr 2 



i?V ) +2— ( / M 2 (j) 2 ) +2 / 



dr 



< / $ 2 (2AKJ) 2 + 2<j\^(f)\ 2 ) (15) 



We select 



61 = -, e 2 = -SUp| Q! | 



and obtain that: 

and simplifies to the following inequality 

Returning back to (|10p , multiplying by 9 and using the harmonic approximation A of A and its initial form Ao 

|AA| + |XA| 2 < \f&\ + 2|XA| 2 + e 
Then for r\ — sup F |Aq| we apply GHI and get that 



i?V ) +2 / tf 2 (^J +2^- / M 2 <j> 2 <2 M 2 cj> 2 



&J.2 



(16) 



(17) 



f|A 



VA n 



and conclude that for suitable 77 > 



if we select /3 2 = C. Also we get that: 



V 2 < J (c 3 (H)rj + 2c 2 (H) V 2 + e) |#?| s 



f d 2 ^\ 2 < [ (\K\+\)V 2 
JF jf 



(18) 



12 /' d ^ 



'J f (\K\+\)$ 2 



"[£) ( 

Furthermore differentiating the equations (|10p we obtain the necessary equations for the higher derivatives of 
the phase function, 

*m = ^, A m = ^-, m = 1,2,3 



dr m dr m 

Specifically we obtain, applying the commutation rules, the following equations 



X*i 



rf 2 *l 

dr 2 



-yA-y$i*i = Si 



X*2 + 
X*3 + 



rf 2 ^2 

dr 2 

d 2 ^ 3 

dr 2 



-yA-y* 2 = s 2 
-yA-y* 3 = s 3 



where the "source terms" are the following: 

si = i? o ^i+R' y'i^+yA i - 

s 2 = 2R / i y*i + —ay, *i +yv 2 

s 3 = - (R'y^ 2 + 2R / 0j0;j y *! + 2R i 0i0 y*y>i +yv 3 -y^) 

Furthermore we obtain the equations for — M = fi£[i + y/i| 2 + K 
T d*i „, ± ^ T ldM 

T dtf 2 T y 1 ^ T lrf 2 M 



(19a) 
(19b) 
(19c) 

(20a) 

(20b) 
(20c) 



*i 



d^ 3 

dr 



•y* 3 = 



ld 3 M 



2 dr 3 



3y*i -y* 2 



(21) 



2R s 0jO y^y s *i - 2R s 0i0 ^y fl *i - r s 0j0 / j >y 



10 



Higher derivatives We derive the estimates for the quantities: 

7j WO = yvi 2 

We obtain, applying the commutation rules, the following equations: 



X7i(0) + 
X72W + 



d 2 li{4>) 

dr 2 
d 2 l2{cj>) 

dr 2 



-yA-y 7l (0) = T 1 
-yA-y 72 (0) = T 2 



where the "source terms" are the following: 

T 1 = - (2R s OOj y0 + R^y 0)y - {r]g o <0^ + t 1 ^^.^^) 

t 2 = -Rm *yv *yv ~yRm *y/> *yv - Rm *yji *y/> + / ^r i y i <t> 

These are the basic equations that will be used in order to derive the neccessary estimates. 
Lemma 5. Let ■& be a cut off along the shell F , then the following holds true: 

c 2 (r,F))V 1 (d(j);F)<V ^(j);F) < C 2 (r,F)V 1 (tf(j);F) 

1.3.2 The shell estimates 

We observe that the equations (fT9|) , ([22|) are of the form as ([19]) and hence are set in the integral form: 

M 2 U 2 ) < / Li} 2 U 2 



(22a) 
(22b) 

(23a) 
(23b) 



l/w-UL 



■d 2 U 2 + — 
/ dr 

where L is an expression comprising /i, R,yR etc. We introduce the quantities 



(24) 



Mr) 
Mr) 



tf 2 U 2 > 



dr 



(25) 



and we integrate (|2"4"]l in the interval I r . £ = ((1 — e)r, (1 + e)r) choosing 

supp(-&)n{r/(r,6)€K}cI r , £ 



and obtain: 



dn (r')\ dn 



dr 



dr 



Ti 



dr' 



l [U;F;§]IL (r')dr' = 
IL (r')r 2 (r')dr' 



(26) 



where 



>*[i7;F;0] = / z? 2 yC/| 2 , 



Ti = / /«? 2 t/ 2 , 
/f 



T 2 = 



Ltf 2 + tf 



d 2 tf 



dr 2 V dr 



dr) 



y?i 2 - ??x^ c^ 2 



We notice first that 



UGALm'+M* 



Ti < C(K) f /" (M) : 



1/2 



D^F;*]) 



1/2 



|/i|(7 : 



1/2 



(27a) 
(27b) 

(27c) 

(28a) 
(28b) 
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In conclusion wc obtain 



I „*) 2 ^« P »/,"» 



Lemma 6. Let -d be a cut off along the shell F, then the following holds true: 

c 2 (r,F))D 1 (^;F) <V°(^;F) < C 2 {r, F)V X (tifr F) 
We will derive an upper bound for tp 2k , tp — logllo then for r G [8^ , 1]: 



1 (.- ,6 1 



'^H |IIol+ w 

Specifically we have by Hardy's inequality 

ip 2kp < C p (2kr(l + e)) p f |V>| (2fe - 1)j V| p < %(2fcr(l + e)f f l^- 1 -*^^ 
x re Jx re 5 s Jx rF 

for < S < 1 that we choose shortly. After Holder inequality we obtain: 

p 

l ^<^l{2kr{l + e)fU l^ 2 *- 1 -^) U |n | 2 V 
Then selecting 6 = fc+e 1 ,_ 1 , e' > we get: 

e p 

4> 2kp <a(j \U A , a k ^(2kr(l + e)r^\C l ,CiV)(k + <'-rr 



Then wc have that 



Summing up we get that: 



^^(.-^(^ 



and hence 



~ /■ W* l^afb ( e'p\(2e'p\^ 

tjx^ (2^)! " [t (2fc!) I fcA*J 

This is majorised after r-function duplication formula by 

oo 

£[C p C(F)(l + £ )r] fe = C (F,r) 

fe=0 

Therefore we obtain that: 

E x 2j 1 x 
(2j)\ ~ 2 C 

I H2<c(F,r) 

Following the usual iteration obtained in the appendix we obtain the usual Harnack inequalities. These are 
culminated in the following 

Lemma 7. The following estimates hold true: 

J 4 2 2 <C 0Q (K)e 

J \^cf>\ 2 6 2 <C i(K)e 

j \^4>\ 2 2 <0, 2 (K)e 

provided that 

' - 2 f lV.il 2 , I \V 2 S\ 2 >e 
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1.4 The lower bound 

The lower bound is obtained by an inductive argument based on the reduction to the boundary of a pixel. 
The one dimensional case indicates the method. Assume that we have a function <f> defined in the interval 
[0,A*]. <M°) = <Km) = ° satisfying 

/■/J rV rV 

c\ I 4>(x) w(x)dx < /i / 4>'(x) w(x)dx < c 2 / 4>{x) w(x)dx 
Jo Jo Jo 

for a positive weight w, < c\ < c 2 . Then we prove min-max principle aalows us to assert that the roots of <f> 
in [0, fj] are at 1 east c\. Then we will derive the inequality through the min-max principle, standard eigenvalue 
and isoperimetric inequalities. The construction is based on the estimates of the preceding section that lead to 
Harnack inequalities for the restriction of the eigenfunction on the boundary of a geodesic pixel. We consider 
first a domain W e C (M. n ,g),n > 2 as is described in the appendix on the harmonic approximation. We drop 
the e subsctript. 

Lemma 8. Let W be a domain with smooth boundary 9W equipped with a smooth metric 7, induced from the 
metric g. Let ip : <9W ->Rfca smooth nonnegative function satisfying the estimate for 7j(V0 = IX^I 2 -' 

c, P°( 7 ^); dW) < V 1 ^^)- dW) < Cjl 2? ( 7 ^); dW), j - 0, 1, 2 (29) 

Let the zero set of ip, N(-0) be (n — 2)-rectifiable. Moreover let (f> : P — > R be such that for r > 0: 

c 30 tV°(4>; W) < V l (cf,;W) < c 31 rV°(^W) (30) 

and 

V°{(j) - ip; dW) < eX>°(V>; dW) 

Then for C 4 = C 4 (c w , Cu, c 20 , c 2 i, c 30 , C31): 

ff°(5W\N(^)) <C 4 

and 

J#> n - 1 (N((f>)) >C vol(W) 

Cq is a numerical constant independent of W, Cy . 

In the appendix we prove that if a smooth function satisfies estimates (|29[) then in every connected component 
of 

W £ = {x e W/(f>(x) > e} 

the following inequalities hold: 

su P < C (F)e, sup iy^| < Ci(F)e, sup IX 2 ^ < C 2 (F)e (31) 

w e w, w e 

for the constants Co , Ci , C 2 explicitly calculated Moreover the harmonic approximation method implies that 
the function n — <t> — d> satisfies the estimate 



|v(c«)| 2 < ct / (wy 

w iw 

The tubular neighbourhood of a nodal set The initial form of the harmonic function <\> denoted by (f>o, 
is of degree m £ y/r: Let 

T e (N(0 o )) = {xe W/\Mx)\ < e} 
and use Hardy's inequality (|34b>|) we obtain 

/ <?< C H e™ / |V0| 2 < c 21 6^r / (C0) 2 < Cre^vol(T £ (^ )) inf H 

JT e (fo) JT,($o) JTe($o) T e (0 o ) 

The usual Moser iteration gives us that 

/ \ 1/2 

sup| K |<Cre™( l —-^ I l«| 2 ] 

T, \vol(T 6 (<fo)) •/T e (? ) / 
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Then according to the usual Harnack estimates we have that if </> c > near y 

T e (0o) = {x G W : \$ (x)\ < e} f]{x G W : 0(x) > r?e} 

we have that: 

sup |« e | < Cre^ne 

T 

we conclude then that <f> ^ <p near N (</>). The harmonic approximation applied on the slice allows us to 
construct the tubular neighbourhoods of nodal sets by the Lojasiewicz inequality for the approximating function. 
Specifically, we have that for suitable choice of e and the tube near the singularities of multiplicity m of the 
nodal set: 

\4>\ < |«| + |0| < Cfae* + l)e 

and selecting re~ = 1, rj = 2 

N(0)cT 3T - m (0) 

The inductive argument For n = 2 we reduce on a disc and then we derive estimates for the zero sets using 
the functions u in order to produce test functions for the application of the min-max theorem, as it was used 
in the Courant nodal domain theorem. We recall here the following lemma from [HSJ 

Lemma 9. There exists rj — rjo(n) £ (0, «] such that with rj £ (0, 770) ifu>i,W2 £ C 2 (B2(0)), \wj\q2 < 1 and if 
\wi — W2IC 1 < 2 then 

U n - l (B 2 _ v (0) f]{ Wl = 0, |Vt»i| > ry}) < (1 + cn)^- 1 (B 2 (0) f){w 2 = 0, |Vw 2 | > |}) 
Estimates on nodal domains and eigenvalues We recall the definition of higher order eigenvalues as 

\k = max mm —7-7 — rr— r 

The min-max principle for the eigenvalue, Afc suggests also that 



A 



k 



mm max — r- — 

s k cm(dw)vGS k yvo^dW)) 

and therefore: 

fT> 1 (v;dW)\ 
Xk < max — — — - — - 
~ ves k \V°(v;dW)J 

Upper bounds on eigenvalues We select as Sk — {(1, ■ ■ ■ , Cfe} where Q : W — > R, Q > defined as follows. 
Let {Ci}^ =1 be the nodal domains, connected components of <9W \ N(-0) and 

C iie = {ieCjMi,NW)>e} 

and 

N j>e (V>) = 5C iie 

We set as 

and 

Olc 4 = 0, j 7^ i, or CilT^nc, = 

and 

\&Q\<r)j 

for Tj,r/j to be selected. After Sard 's lemma Nj^ip) is smooth for suitable e > and hence we assume that 
dCi is also smooth We approximate harmonically tp and construct a smooth partition of unity, each member 
supported in a connected component. The Harnack inequalities in the appendix suggest that 

\$-ip\\2,9w<C(dW)e 
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and 

sup \ip\ < r\ 
aw 

We compute that 

A fc < E -^ 



where; 

aj = vol(N j;e (^)), bj = vol(C Ji£ ), bj = vol(Cj,e \ C J>£ ) 

Sard's lemma again allows us to choose e, /it so that 

bj ~ M&j 
Furthermore the isoperimetric inequality suggests: 

.. n-2 / 

6™- 1 < C Uj(l + ej + sup |yV>|) + bj sup |/i| 

V 8W 9W 

Moreover for suitable e > 0: 

fc 
^ 6j = vol(W) - vol(T e (») > (1 - e)vol(W) - (vol(N(»))^ 
j=i 

Let 61, 6fe be such that: 

61 > -vol(W) 
k 

and similarly 

6 fc < ivol(W) 

We conclude 

fc(jr ri - 2 (N(^)))^ 



A fe < 



vol(W) 
Number of nodal domains The first eigenvalue of the laplacian in C^ e satisfies : 

Al(c ^ } - j^v 2 

Therefore 

Al(Cj, e ) < C01 

We need now a lower bound for the first eigenvalue of Ci ]£ : we select e, so that we avoid dumbbell shape of 
the nodal domain. Checger estimate of the first eigenvalue combined with the isoperimetric inequality on the 
spherical piece suggest that 

2c i > Ai(Ci, e ) > — - inf \h\ 

vol(C,, £ )— c s „ 



vol(Ci l£ ) > IcoiH- mf \h\ 

Hence since at least one nodal domain should have volume at most 

c^ < ^vol(F) 

hence we have that 

k < vo1(F)cq7 1 
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Lower bounds on eigenvalues Wc split <9W in k domains, through the use of harmonic polynomial ipo of 
degree d(k) ~ — Jj — - as follows. Let {r]j}°^ be a decreasing sequence of positive numebrs and 



and 



Then wc exhaust 



where 



T„(N$o)) = {x e 8W/\Mx)\ < rj} 

w j = (t„,(n$ ))\t w+1 (n$o))) f|aw 

DO 

F = W_i (J Wj 

3=0 



2 



W_i = {x e dW/Mx) > Vo} 
Then we have for a smooth function <j) '■ W — *• R and j > 1 and for N > 0: 

The definition of Afc along with the preceding inequality gives that 



Select the sequence {rjj} so that 



Therefore we have that 



Finally we conclude that 



vv / \— " ° Jw 3 



00 f 2 ^ f 

E / ^ 2 <^E/ <? 



r/l'A-i 

Afc > 5 



1 + V * W 



Modification of the pixel for the singularities of Hopf 's strong maximum principle for <\> guarantees 
that N(^>) meets transversely the boundary of the pixel. The comparison of nodal sets reduces the problem to 
the estimation of the nodal set of (f>. Therefore using the geodesic spheres starting near N(-0) we obtain by the 
coarea formula, away from singularities of <j>: 

jr™- 1 ^ (c/>)) = E / dr / d.Jf n - 2 (N(iP)) > ct- 1 ^ 

The eigenfunctions The eigenfunctions fullfill the hypothesis of the preceding theorem and therefore 
Theorem 10. Let u : (M n ,g) -> R, A g u = -Aw,P C (M n ,g) be a pixel. Then 

jT'-^n^a) n p) > C(P)S\ 

We are placed in the P with boundary spanned by the fronts F^, £ = 1, . . . , rrii 1 ...i k : 

n+l m n-i k 

^U U U ** 

fc=i allu,...,j fe fcl 

The boundary of the pixel is not smooth, but we apply the smoothing method described in the appendix. Near 
each front we have the representation of the cigenfunction in the form: 

u(r,6) = e-7S" (r,fi)+/9aa(fi) siii(^i0(r,e)) 
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In the preceding paragrpah we obtained the estimates that <fi satisfeis with the constant en ~ A. Therefore we 
have that inside a pixel of diameter -k= : 

Jt? n - l (N(ux) n P) > C(P)X-^ 
and since the manifold is splitted in A 2 pixels of trhis size we have the required estimate. 

1.5 The upper bound 

[D] proved that the Hausdorff measure of the nodal set contained in a pixel P with its boundary smoothed out 
is majorised by: 

J^fNWnP) <- f |Vlogg| + \/^Avol(P)+vol(<9P) (32) 

2 Jp 

where 

q(u) = \Vu\ 2 + -u 2 
n 

We split the set P in three parts 

• The part that is free of nodes: 

P e = {x e P/\u(x)\ > v } = P \ T„(N(u)) 

• The tubular neighbourhood of the nodal set T v (N(u) n P) is splitted further as 

T„(N(«) n P) = #,(N(u) n P) u y n (N(u) n p) 

The regular part of the nodal set 

^(N(u) nP)={ie P/Hx)\ < 77, \Vu(x)\ > t]} 

and the neighbourhood of the singular set: 

Sing^(N(u) nP) = {ie P/\u(x)\ < V, |Vu(jc)| < 7]} 

The problems in (|3"2"j) arise in the singular part of the nodal set. We will estimate the behaviour of |Vlng| 
near Sing e (P). We will use induction with respect to the multiplicity of the nodal set, introducing new pixels 
of multiplicity bounded from below and approximate harmonically the eignefunction there: 

m(\) a 

T,(N(«)nP)= |J |J Sing Am 

where Sing^ m is a connected component of mulitiplicity £ 

Sing^ m = {x e V/u{xf + ■■■ + |V'-y:c)| < V, |V«(a:)| > v} 

Notice that 

Sing , m = P ei Sing lm - Reg„(N(u) n P) 

The harmonic approximation of u\ in Sing^ m is denoted U£ >m . We introduce the localization functions Ce,m 
with supports supp(&, m ) C Sing^ m We split the integral as 

\Vlogq\<J2 Co,m|Vlogg|+ J2 / 0,m|Vlogg| 
Therefore we set 



hACm] = / OjlVloggl < C f CijV 
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where we set _ _ _ 

IVIVmI 2 ! + c|u||Vu 






|Vu| 2 + cu 2 ' > 3 7 S 



We will estimate the terms appearing in the right hand side of the inequality. If £ = we have that for 
D(r], A) = r f it \ qn+3 where t is a parameter to be selected: 

f] < |w(x)| < D(r],X)rj 
|Vu(x)| < D(t],\)V\t) 
|V 2 u(V)| <CD(t},\)\ti 

Therefore we have that 

IUoj < A-* +9n+ V 
and we select as 

r/ = X , m = — I 

s ns 

If £ = 1 then we split the tube around the regular if of NV(u)) in the layers 

DC 

•^ = U ®i' ®i = fa e /^ - l«( a: )l^~ 1 '?' l Vu l ^ ^ 
.7=1 

We have that in ^ 

|u(a:)| < D{rt,X)B j r) 

r)< \Vu(x)\ < D(r],X)VXr) 
|V 2 u(a;)| <CD(r},\)\ri 

In this case we get 

IU 0J < X-S+vn+irT' 

and we select 77 as before. If £ > 2: For this we need sharper estimates and we appeal to the harmonic 
approximation in order to use Lojasiewicz inequalities. We exhaust Sing^ m : 



Sing, >m = |J Qi(6, r?), Qi(9, r,) = T{ m \ T^ 1 



3=0 

for 

TJ iOT = Sing £m p|{x € P/^' +1 »y < |u| < 0*17, ^' +1 »y < |V«| < Fr?} 

We introduce in each Q £ m (0, 77) then 

and drop the indices 

V|Vu| 2 = V(Vk+ • Vk_) + V|Vw| 2 = (V 2 k+) • Vk_ + (V 2 k_) • Vk+ + V|Vu| 2 

The harmonic approximation method combined with Harnack inequality suggests that in Q £ (6 , 77) , £ > 1 we 
have Bernstein's inequalities 

|Vw| < D(r], X)\fXr)d j |V 2 m| < D(r], X)Xr]d i 

and also 

|V/s_| < CiA0 2(j+1) r7 

Set 

/3 = D(ri,X)r]6 j 

Furthermore since \a + b\ > (1 — e) (a 2 — ^6 2 ) then : 



«r >(i-e)«- -M > 1 - ^^ « 
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and 

|Vu| 2 > (I - .-) I 1 - — ) \Yu\ 2 > i 1 - :) j I - — j \u 



-' > n - -. ( 1 _ ^M I IV«|2 > n -,W 1 - ^M 1 |ft|2(i-^)) 



In Q^ m we have: 

u 2 > ( 1 - ^^ ] 2 W+V , |V«| 2 > (1 - e 2 ) [ I - ^^ ] ^(i-, W )0+i)^2(i-, 



V ^ 1 n2(i+D^2 irr..i2 ^ /i _2i / i V A/3 \ /1 2fl-i/mifi + 11_2fl-v(f)) 



The integrand is estimated according to the estimates derived above: 

We have that due to the multiplicity bound we ahve that 

vol(QiJ < CX-^VXrjP, D( v , A) < C(n s \ q ) n 
We set rj = \~ m and compute 



Aqn+ 1 

+ 3 



Summing up we get 



and therefore we have that 



2(2sn + 2iv+ 1) 

v OO ~ 



cvx x _ 3 



3=0 j=0 Js e,i 

cVa 



e 2 »' 



We will use repeatedly the following method that we call harmonic approximation method. The domains 
encountered here W have boundaries with normal crossings singularities: the singular set J^(dW) is given by 
the transversal intersection of hypersurfaces: geodesic spheres with local equations s\, . . . , sg,£ = n, n + 1. The 
piece of the hypersurface Mi = {x £ W(si ■ ■ ■ si + e)(x — 0} near S* is for suitable e a smooth hypersurface 
close to ^(ffW). We will consider the domain W obtained by replacing the singular part y(W) by Jf e with 
repacing the defining function through cut-offs by the function given there. Let F : W — > R be the solution of 
the boundary value problem: 

AF = 0, F\ d ^ = F 

Harmonic polynomials We will also approximate the harmonic function defined in the pixel F by a sequence 
{-FiijnGN of functions such that 

Ao^ = 0, A F„ = - ]T TV* f^pi - g^d^djF^ (33) 

C/ <J-"i \JJb •) 

where 

9ij = % + ^»j) = diam(W), if) = - log (3), ,9 = det(&j) 

and for j = 0, 1,2: 

||v j '-k|| <Chq 2 - j , 

Integration by parts after multiplication by C, 2 F n and incorporation of the preceding estimates along with 
Young's inequality leads to: 



/ C 2 |VF„| 2 < o 2 / C 2 |Vf„_ 1 | 2 + c 2 / (|vc| 2 + C 2 K 2 

Jw iw Jw 

supp(|VC|) C A(W) = {xe W/d(x,dW) < e} 



/w Jw Jw 

and 



and 

IV'CI < § 
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We select C\ip 2 — 1 then 



Similarly we have the inequalities: 



/ C 2 |V^„| 2 <C I |V(CF )| 



and 



/ C 2 |V 2 F„| 2 <cfp 4 / aV^-xlW / (lACI + lVCI^lVF^I 2 ) 
Jw V Jw Jw / 

/ C 2 |V 3 ^| 2 <C 2 f P 4 / ( 2 \V 3 F ni \ 2 + P 2 f C 2 |V 2 J^_i| 2 + p 4 / (|A c | + |VC| 2 )|V 2 F„_i| 
Jw V Jw Jw Jw 



Therefore we have that after iteration: 



( 2 \V 2 F n \ 2 <C C 2 |V 2 .F | 2 
w Jw 



and as well as 

/ C 2 |V 3 F„| 2 <c/ C 2 |V 3 F | 2 
Jw Jw 

The Nash-Moser iteration that we describe in the sequel allows us to bound the sequence in C 2 (W). Rellich 
lemma allows us to extract a sequence that converges in i? 1 (W) and therefore we use it in order to approximate 
the initial function by a harmonic polynomial with any accuracy we desire. 

1.5.1 The brick localization details 

Let P = p(^;*i'---' l fc) be a brick of size determined by the parameters (r(P), /i(F/). We use coordinates x — £ + c 
for c denoting the centre of the brick. We construct cut-offs Let ip <G Cq°(R+) be the following function: 



Mt) 



1 if o < t < 1 

if t > | 



Then the following function that localizes in the brick: 



U 



<H*m 



Also we will use the function that localizes in the neighbourhood of the zeroes of the function 

h:V ^R, 

N(ft) ={ieP: K(z) = 0} 



then the function 

localizes in the set: 

We can prove inductively that: 

and 



t(X) = to(xo)ip 



h(x) 



N e (h) = {x<=P/\h\<e} 
|V^£o| < Cj 



|V J '^(x)| < Cj J2 \V ll h\ Pl . . . \V im h\ Pm \V jl ip\ gi . . . \Vi'iP\ qi 

where we sum over all indices iipi + • • • + i m p m + jiQi + ■ • • + jiQi = £ and hence that again for indices 
HPi H himPm =£,£ = 0,...j: 



|V^(£)|<Q^ 

These will be used succesively in the sequel. 



h 


pi 


h 
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1.5.2 Lojasiewicz, Hardy for functions of the form /ioJVj 1 

Let hg be a polynomial function in rectangular coordinates in B l R then we have the following immediate result 
Nodal volume 11. The function h = ho ° N^ 1 is a function that satisfies the following 

• The mulitplicity strata E m (/i) of the variety N(/i) are mapped to $(E m (/io)) = E(/i). 

• T/ie Lojaziewicz inequalties hold true 

\Vh\ > CilKl 1- ' 1 , |fc(»)| > d(x,N(h))) 
The first conclusion comes form the chain rule in many variables: 

D a h = S^ )7( ,- ) (i^)(*( i r))(iy»*) ei • • • (IT"*)* 
where thje sum extends over all multiindices a, ft, j(j) G N™, j = 1, . . . ,£, e\, . . . , e# € N such that: 

\/3\ = l,e 1 - fl + --- + eae = \a\-\0\ 

Exchanging the role of h, ho, $, N we get the defining equations of the equimultiple locus. For the second we 
compute: 

|v/| = | J D$,(vK )($(x)| > c\vh ($(x)\ > c'\h^{x)\ 1 -^ = c'lf&i 1 -^ 

where we have chosen R, r so that 

\D<f>(x)\ = \A + B(x)\>\A\-\B(x)\>^ 

where \B\ < ^. Similiarly. 

\h(x)\ = \h ($(x)\>c 2 d($(x),N to Qy* 

now due to the first inequality and the definition of the exponential map we conclude that: 

d(*(as),N(£o)) > c'd(x,N(h))) 
A consequence of this is that (GHI)'s hold for such functions. 

1.5.3 Hardy's inequalities 

Let P : R™ -^Rbca homogeneous polynomial of degree m and N (P) its set of zeroes 

N(P) ={ie R n /P(x) =0} 

Moreover let / e C °°(R n \ N(P)) there exist constants, [P^ < c 3 {H) = t-^w + 0(e)n > 2 



(n-2) 2 



/ wr-f^c^H) j iv/i 2 

f < c 2 (H) f |V/| 2 , 
JR" 

f < c 3 (H) f |V/| 2 , 

JR" 



P 

AP 



(34a) 
(34b) 
(34c) 



From the euclidean Hardy' s inequalities we obtain the riemannian versions by modifying suitably the constants 
by 1 + e. 



21 



1.5.4 Integration formulas 

Let T be a tensor field of typr (p + 1,0) then we introduce: 

-^■l.-* )ii...i p k — J- ii...i p \k ^-i\...k\i p 

-D(r)u...» P -i = g 3 T^...i p -ii;i 

The general integration by parts formula reads as 

Lemma 12. Let T be a (p,0) tensor field on the riemannian manifold (M™, g) and <p — \T\ x, x> a smooth 
cut-off function supported in the domain K. Then we have that 

f 2 |VT| 2 < I l -\A(T)\ 2 + \D(T)\ 2 + I 4? ^(Rm *T*T) t + f \W X \ 2 \T\ 2 

We set T\j = T^..^ _ x j and hence we have that 

4> 2 \A{T)\ 2 = 2(f> 2 (|VT| 2 - T- lk .jT' lj ' k ) 

The last term gives that 

4> 2 Ti k ;jT^ k = (^TikT^).. - (f?T lk TVM - 20T ife ^T^ fe = (BT)+(I)+(II) 

Furthcmore we have that 

v 
(1) = (j) 2 T ik T' lj *' J = </> 2 T ifc (D(T)i). k + cj} 2 Ric k s T is T ik + (f> 2 J2 Rm l s ekj T h T ik 

for \i = i\ ■ i£-isi£ + i ■ i p . The first term then is written as 

<f> 2 T ik (£>(T),). fc = (ST) - 2 |£(T)| 2 - 2m^j(D(T)y 
The term (II) is written using that (j> = x|T| fc_1 : 

(II) = 2<\>T lk <\>,A(Tf k + iviogx • V|T| 2 + (k - 1)^ 2 |V|T|| 2 
In summary we have that: 

jV|VT| 2 < J^ 2 (\D{T)\ 2 + l -\A(T)\ 2 ^ +cj^ (|Rm| + |R*c| + |Vlog X | 2 ) cf> 2 \T\ 2 

1.5.5 Functions 

For a function / : K — > R, supp</> C K we have that: 

2 |V 2 /| 2 < C f cf, 2 (|A/| 2 + |Rzc||V/| 2 ) (35a) 

K JK 



I 2 |V 3 /| 2 < C [ <f> 2 (|V(A.f)| 2 + |Rm||V/| 2 ) 

JK JK 



1.5.6 Curvature 

The diffrential Bianchi identities are written as: 

,c«ff(Rm)i,fci; m = Rmij m i- k 

also: 

&tf(Rm)ijki;i = A(Ric)j k i + A(Ric) ifc j 
We recall here the Bach tensor: 

Byfe = A(Ric)ijk + , —ry, _ . [gijVkR - ffifcVjR] 

and the contracted identities are writen as 

M{Ric) l =y^R 
Therefore we have that 



5 

I 2 |VRm| 2 = 4 I \B\ 2 + -? ^|VR| 2 + ^(Rm * Rm), 

JK Jk (n-ljtn-zj 



(35b) 
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The iterative method Let g, \ '■ W — > R be smooth functions and fc:K->Ra polynomial weight function 
of degree m then set 

Wj = {x e W/ 6>(1 - 6> J )| < |K(»)| < (1 - + 6> j M 

and suppxj C Wj . This for instance is given for 



Xj{x) = e\jY 



h(x) 



9(1 - 6 j )ri 



+ di)vl \ h{x) 



We suppose that the smooth function g satisfies the inequality, for positive constants 7 > l,e — 2,4 and any 
smooth cut-off x : 

x 2 |v ff | 2 < 7 / x 2 lffl e 



w 



w 



Then Sobolev inequality suggests for s — ^zz, 1 < P <2,k,£,d> 1: 



p/s 



w, 



d\u\-£„k 



X d \h\ 



<C Q k (\h\- e x d g k - 1 ) \Vg\ p + i 



V/i 



/1 



(x d N-y) P + (36) 



+p / x {d - 1)p \Vx\ p 9 kp \h\- pe (37) 

JWy-i 



The first term then gives for r = tt^- 



(\h\- e x d g k - 1 ) \Vg\ p < I / 



cip 



(ffi-y- 1 



p/2 



W.,-i 



/ p |v 3 | 2 



< 



< C x ((l - + ^')? 7 )^7 f vol(W J ) (*-i")-° / X ^l^r rf 3 (fc " 1)r 



w, 



(38) 



since we have 



p/2 



;.><.;• 



w, 



X 2dp \Vg\ 2 ) < 7? ((1 - 6 + 0*» "- 1 vol(Wj) (*-d- / X^ £ |^r r£ 3 (fe_1)r 



Wj-i 



The middle term after application of (|34bl) gives that: 
p 



w, 



Vh 



(x d \h\- e g k ) 



w, 



x*^*- 1 ^ 



w, 



Vh 



X dp 9 2 



< Ci((l - + ^"V vol(Wj) ( %-V / x^|/Jr rf g (fc - 1)r 



In an analogous way the last term leads to: 



\^x\ p (h' i g) k < Ci((i - e + e i )7 1 y l ^vo\{w J ) %-iV 



w, 



X^\h\- ri 9^ 



w,- 



In summary we arrive at 



(x d \h\~ e g k ) ] <(7 3 /c*(?7(l-0 + J ))^7*vol(W J ) i T^fe L 

vv, v ' I \JWj_i 



(fc-i+f ) 1 

fe 'r(fc-l) 



x 5<lr|^|-r^(fc-l)r 



Notice that 



r = qs, q{n,p) 



1-2- 

n_ 

1 2 



2 2 

1 < ?(ra,p) < 2 

n n 
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We conclude with the basic inequality that we will iterate 



where 



— I (~iks 



<C* 



1 



G (k- 



s n 

r = -, a= -, 

a n — 1 



Vj-l JWj_i J 

Vj =vol(W j ) 



Cj — C3 



k V sh ^ lt {e{i-e j - 1 ))' 



k fc-i 



<-y 2 k y . 



ft 



r(fc,-l) 



"j+1 



and we replace j3j by the upper bound 



„a(„-2) 
ft < ^^ = <* 



In order to bring this to the standard iteration form we do the following: 

H.J + 1 



kj — 



1 



Finally we arrive at the basic iteration inequality: 

1 



Ij+i - 



G a 



<d 



y v j+ i j w . +1 
Then we have the iteration inequality 

The iteration leads to the inequality 



G a 



U J JWj 



I ] + l < Cjlj 



SU P \ G \^ D \ uw ^ / GP 
Woo V vol ( w o) 7wo 



D= lim (Ci...Cj) 



(39) 



We select the local density parameter as 9 ~ -^ , t > 0. The constant is estimated through elementary inequalities 
of the form: 



7'a 2 



1 ~ l / 1 



We arrive finally at 



3=0 



c = rj 



I 



ytj J /yt _ ]_ jt a 2 _ I 



fnn 4- a+1 P"(t + 1) , t-1 
' t+ o(n-2) ry 21 ^^^2 



We will denote the constant in the form: 



D(V, 7) = fa'7T, s = ^P + 



o + l 



9 = 



p(t+l) t 



na(n - 2) ' '' 2i ' n(n - 2) 

Back to harmonic approximation The harmonic approximation method suggests that 

sup\ K \<D(r 1 ,\)(— i— /" 

w^ V vol ( w o) iwo 



1/2 



Similarly we have the higher order inequalitcs 



f 



su P |V«|<£(77,A 2 ) [l + D(n,n)X- . 

Woo V von Wo J iw» 



|Ric| 



1/2 



VOl(Wo) 7wo 



1/2 
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Let u > be a harmonic function. Then set for j = 0, 1 



hj = 



V J+1 i 



\\73u\ 



Hn = 



VlVul 2 



iv« 



and compute 



or that 



A g ho > -(\Ric\ + 2)hl + -hih > -(\Ric\ + 2)h 2 



f C 2 |Vfe | 2 < / (|R* C |+2)C 2 ^ 
Jw Jw 

We note that the integral is of the form 

/ <fhl<{l + e 2 )c 2 (H) f <f> 2 \Vh \ 2 + \\Vcj>f 

majoriscd after application of Hardy's inequality. We compute that 

1 



|vft | 2 <(i + £ )W + -M 



then we find 



</> 2 |V/i | 2 <(l + e) 2 / cj> 2 (H h 2 + ^h )<(l + e 2 )c 2 (H) f |V0| 2 ( \ + h 2 
w iw e Jw \ e 



If n > 1 , e' > then select 



n 2 — 2n — e' 2 9 n 2 — 2n 

e 2 < 



n 2 - 2n + 8 n 2 - 2n + : 



and conclude that 



and hence 



w 



^iVftol 2 < c„(e')y w |V0| 2 [ 72 +h 2 oj , c n (e') = - a ^ -g 



Repeating the procedure we get 



2 ^<c„(e') / ^(l + e'/i 2 ), 0i = |V0| 
w iw 



2 /io<c„(e') / (l + e 2 /i )W 
w Jw 



and 



4ho<c n (e') / (l + e 2 )|0 3 | 2 , 03 = |V0 2 | 
w Jw 



We apply this formula for <f> = (|Ric| + 2) 1 / 2 £ that 



and we obtain 



where |V£| J ' < 07 J 



3 

/ C 2 |V^ | 2 < / E c„(e') J |V^Rz C ||V^| 2 
Jw Jw i=1 ^ . 

suph <D(t,-f- 3 ) [ |R|, t = max(ti(W),t 2 (W),t 3 (W) 

w Ja(w) 



^IV/zoI 2 
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Growth of a function near its nodal set We assume that 



w iw 



2„,2 



Let 

and 



W,- = {xe P/(l -6 j )6r]< \u(x)\ < (1 -0 + e J )f]}C\B c ^_ 
A(W,) = Wj \ {x G P/(l + ^ +1 )fl 2 < \u(x)\ < 9(1 - P)T,}r\B Ct j_ 



The cut-off function £ satisfies the following estimate: 



|v £ CI < % 



We apply Hardy's inequality 



CV=/ w---(Cw) 2 < Cl (ff) sup|u| / |v(Cm)| 2 < 

w, Jw 3 \ w ^ / j Mj) 



< C sup |u| (1 + e) 



w,- 



-/ |VC|V+r/ «) 2 



Therefore we have that close to T,, (N(u)), for 77 = (2ci(iJ)r(l + e)) 2 . We have that 



cv < -i 



|VC|V < Cl 



r-e6» 2 7 A 



(W, 



We have that: 



u 2 < 2 / u 2 + k 2 

A(Wj) 7a(w 3 ) 



and the second integral is estimated again as 



A(Wj) JA(W 3 ) 



and hence 



U 2 < 



'a(Wj) 2 — 9» Ja(Wj) 

The coarea formula suggests then after Lojasiewicz inequality that 



c 2 « 2 = 



A(W ; ) 



(e-e^ii+Ov 



dp. 



(6-6i)r) 



, f (i-6+ei)v 



{u=n} 



Nu 



dfi 



(l-6+8i)£ri J{u=p,} 



Vu 



< 



/i u+1 a(/i)dii + / n" + a(fj,)d/J, 

{6-63)7) J(l-e+$i)£ri 

Inside a ball of radius r applying Crofton formula if the multiplicity of u is m 



a{n) 



{u=fj,} 



da^ < ct 2 m 



Therefore we find that for ct,- = 6 — 6 J : 



J C 2 . 2 < t-^V+ 2 f^i^ , (i-^a-^) 



'A(Wj 



1 - 2ct, 



Finally we have that 



C 2 u 2 < Ct — ~r] 



v+2 



W, 
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Moreover we recall the folowing identity from [P]: 



v v^-v / C v 

art Jw 



for Q = |Vw| 2 > c\u\ 2{l ~"\ This leads to 



rfri—ir* / C ^ < 

G"7 Jw Jw 



VQ 



w 



C a u"+ 2 < 



VQ-Vu 2 ^ 



w 



VQ 



Q 



C 2 u v+2 



We apply Hardy's inequality and get 



for 



Finally we get that for r\ > t\q\ 



and 



r)I'(v) <c(l + r _2 *V) tv v I(v) 



m = »r 8 / C 2 u 2 

Jw 



I(jj) < Ce 

( z u 2 < Ce 



C„( l+T 2 n" +L )TT] 



I(Vo) 



w 



, Vo J J w 
Morrey estimates Let e < 1, 0<7<lor7<0, p < 2: 

u e = y u 2 + e 2 , ij) e = \ogu e , w = u\ 
and for (, supp(C) C W: 

2 2 



C 2 u 2 



f C 2 |Vu e | 2 <r/ C 2 u 
Jw Jw 



Then for q = -: 



I 2 ,,,? 



/ \Vw\ p C <Cx{t) I |VC| 2 ™ 
Jw Jw 

I |VV> e | 2 C 2 < C 2 (r) / |VC| 2 + C 
Jw Jw 

The inequality (|41a[) follows after selection for £ as uj _1 gives since 

|Viu| =7w 1+ t|Vw £ | 
that 



/ 

JW 



C 2 |Vw;| 2 < Co ' 2 - 



7 2 Jw 



The inequality (|41bl) requires the additional assumption for r > 0: 

C 2 |v 2 u| 2 <r 2 / cV 

/w iw 

We start selecting values i>i , . . . , v m > and assume that 

u = Vj + hj 

making the following choice: 

\hj\ < e\vj\ 



(40) 

(41a) 
(41b) 
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then 



fe+/b) 2 >(i-e)k 2 -- >&-#% 



/; 2 
e 



A2.2 



We approximate harmonically hj in suitable bricks selected so that we use the initial form of hj . Hence we 



have that for ip = log(u e ),-0 = log(/i) 



p jp 



2^2 



Therefore let h be the harmonic approximation of h in W and we set: 

h = h + k 
The standard harmonic approximation method estimates combined with partial integration leads us to 



/ C 2 |V 2 k| < D( V ,r 2 ) [ ( 2 u 



The estimate of the preceding paragraph 



|Vk| < sup|V«| < c / (V 
Wo Vw 



1/2 



D(r 1 ,T) + D(r ] , f i)T\\Ric\\i 
vol(WO) 



We compute for e < 1: 



uf=u 2 + 2ku + k 2 + e 2 > (1 - e 2 )u 2 + (1 - — )k 2 + e 2 



Then we select e such that 



and therefore we get that 



^l-^) + e 2 >\e 2 



N< 



^(l-e 2 ) 
Let then m denote the highest multiplicity of u. We apply the preceding estimates to conclude that for 



2(1-^) • 



c : 



2|VW £ 



< 



1 



w 



2(1 - e' 2 ) 7w 



v«, 



«;', 



^'+'"/ C>r I $r 



2„,2 



Now we use the estimate of the preceding paragraph 
and obtain that 



C 2 |V^ e | 2 <C e / |VC| 2 +ct-^ v »+ 2 I Xr 



w 



w 



w 



We set efe(m) = 1 — — and then 



C f (u-^o 



w u Jw u 



eif\2 



<Cl(H) 



VS 



The constant in Hardy's inequality is ci^{H) 



\V^\ z u- 2 ^C<{l + e) Cl {H) 



(d-l)^' 
1 



W 



d>3: 



(«--C) +/ ^ 2C1 |VCI 2 

iw 



w 



We select then e so that : 



Hence we have 



«- 2ei |VC| 2 +e 2 / |VVfw- 2ei C" 
6 iw iw 



e< 



(d - 3)m + 2 
(d + l)m-2 



C 2 |v^| 2 < (1 + £)ci(g) 



w 



u- 2e \V(\ 2 



w 
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and we conclude that: 



f C 2| V?/ ,|2 < {l+t)ci(H) J s - 2e|vc|2 (42) 

Near the zeros of u, u, i.e. selecting r\ suitably small then we retrieve the same inequality. Further treatment is 
required when we introduce in (|42p 

and we get: 

' W|vvf < ci 2 f u^ (Vv ( *- 2) |vvf +V 2( ^ X) |W| 2 ) 



/w ./w 

through the elementary inequality for x € (6* , 1]: 



2 / , ,S , _2 4~^ _1 



|logx|< - I \x\v + \x\ ^ < s j 

We obtain for g < 1, <5 = (e^) _1 : 

tfV|VV| 2 < 2C£ 2(e ^ +1) / ^"C 1 -*)) [tf 2 |V^| 2 + V 2ee£ |W| 2 ] (43) 



w Jw 



The iteration for the lower bound. We follow the method of |SL2j . We select as £ 

v 2(i-l) § 2(al-b)^ v = ^ e - a , a>l + b 

obtaining that: 

/ ^2(Q£-f,) v 2(f-l)| Vu |2 <^2(e e £+l) f y 2ti-A | ^^2(^-6) + ^(£-1)^2(^-6-1) 

iw iw 

We get that 

/ |V(| V |V^ b )| 2 <^ (ee+3) V0l(W)+^ ee+1 ) f w 2^2(a£-fc-l) 

Therefore we have that for a — kb, b = , m , m = -J-* 

1 fe— m ' £—2 

(^)/„ (l " |i,)2 ") 4 <^<««>+^«> (^/ w (l*") 2 ' 

Hence we have that 
Selecting a sequence {£j}: 



tj = ^' ii = ( ~ I (\v\v) 2kJ 

Hence we obtain 

Ij+i < k Jiee+3) + k j{ee+1) Ij 

or that 

Ij+i <fc (ee+3W (l+/i) 

Examining separately the two cases: for some jo : 

ho < 1 

and the complementary case we arrive at the conclusion 

ij < (2k j °y- jo 



2<> 



i/fc J 



Therefore following the reasoning in |SL2j we conclude that 

f e cihk-ff| < c*vol(W) 

which implies that 



Wo 



(u 2 + e 2 ) + / (u 2 + <?)-*■ < C 2 vol(W) 



and for small p > 0: 



w 



inf u e > C I / ii' 

Woo V^oo Jw, 



We follow again |SL2] appealing to (JUT) and get the bound: 

inf Id > C | — 



W„ 



'oa JW 



1/P 



1/p 



for anyp < ^. 

1.6 The two dimensional case 

We will derive a version of Hardy's inequality for the two dimensional situation that is not covered in the general 
case. Therefore we start with the radial blow-up of the plane covered in the 



We set in C\ 



d = R 2 \ {|xi| > e\x 2 \},C 2 = R 2 \ {In | > e\x 2 \} 



r£ r 

X\ = , = , X 2 



and interchange the roles of x\,x 2 in C 2 We obtain for 



P(x) = r m R^), J = 1,2 



the elementary identity: 



| V p|2 = r 2(r„-l) 



^+(i+av|((i+a-^i) ! 



Therefore we have that: 



Similarly: 



R 2 



< 



R 2 



P2(l-i) 

|VP| 2 ~ m 2 + {{1 + e)R' - m£Rf ~ m 2 
VP 



2 2 

771 



». i? 



(l + ^)__ m £ 



We compute then that for 5 < 1 and / G C£°(R 2 \ ({P = 0} U{log \P\ = 6})) then 

/ I f 2 <^L[ IV/I 2 

y R2 p^ii og m| 7 - uogwi y R J Jl 

and 



IVPI- 



;/ 2 < 



C 2 (P) 



/w P 2 |log^| " |log(<5)l Vr2 
Then we start localizing in areas of constant sign for R: 



IV/I 1 



h* = 



R 2 |fog^! 



VP 



p 



XL/' < Ce 



x,/ 2 



R2 r|logfr JJ | log 



(! + £')— -m£] Xif 



We used the inequality: 



(a + fe) 2 >(l- e ) 2 (a 2 -^& 2 
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alternatively in the regions (logi?) 2 > ?7 2 (logr) 2 and (logi?) 2 < 2?7 2 (logr) 2 . The inequality is majorised after 
integration by parts in the radial variable through the elementary inequality: 

1 2 , 4 f 00 2 



/ r 2 |logir -|log<5|7o 
This proved easily by splitting the integral after arranging e according to the support of g 

rl + e 

1 /■ 1 f 00 1 

g 2rfr = / "m T\9 2 dr + / — —g dr 



, r 2 |log|r 7o r^llogjr ./,- (1+ .> r 2 |log|| 

The first integral is written after integration by parts: 

„ f Sl+ ' 1 , 1 2 1 /*°° /„ 1 + e f 5l+ ' 1 ,2 

2 / "T - T^ + -77, J77,9 < - / .9 2 + pj— T7 / — ^/ 2 

7 rlogfr r 2 (logf) 2 e J | log<5| J r 2 1 log ^ | 

then we choose s = [ °| [ and we get: 

21 T .9 2 < -Tj— FT / 5'2 

/O Hlogdr e| log J| J 

Similiarly for the other integral, setting as g 2 — rf 2 and splitting the integral in two pieces then 

1 ., 2C r f°° .„ 2C f f°° 1 



rf<j^ rf' 2 + T ^- —±—rf 2 < 



r 2 |log|| J ~ | log S\ J J | log 6\ y r 2 |log(f)| 

The two dimcnional inequality is obtained by a direct application of the usual one dimensional inequality in the 
£- variable after the formula: 

B!\ 2 J„ ^ m 2 



n) * a l c+ £«-*> 



■j. 



1.7 Curvature estimates 

In the integration by parts formulas we substitute the curvature identities we find that: 

/ 2 |VRm| 2 < 3 / </> 2 (|R| 3 + |Rm| 2 ), f 2 |VRic| 2 < 3 / 2 |R| 3 (44) 

jk jk jk jk 

The iteration scheme suggests that 

sup |Rm| < D(r], fj,) inf |R| 

w * w* 

sup |Ric| < D(t], h) inf |R| 
w* w * 

1.8 Local properties of eigenfunctions 

In the case of an eigenfunction we have the following 

<f> 2 \V 2 u x \ 2 <\(\ + D( V ,n)\\Ric\\ hW )([ cj> 2 u\\ (45a) 

2 |V 3 u A | 2 < A 2 (A + U(7 ? ,//)||Ric||i, w ) / 4> 2 u 2 x (45b) 

K JK 

Performing partial integration to the term: 

/ 2 |R JC ||V.f| 2 = -2 / |Ric|/W-V/- / 2 /V|Ric| ■ V/ - / |Ric|^ 2 /A s / 

JK JK JK JK 
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Young's inequality along with harmonic approximation for -\/|Ric| 2 + e leads to 



<j) 2 \Ric\\Vu x \ 2 <CX \Ric\(j> 2 u\ 

K JK 



Similarly we get for Finally we conclude that 



^ 2 |V 2 WA | 2 <A 2 / ^(l + ^\ (46a) 



6 2 |V 3 UA | 2 <CA 3 [ 02 (1 + 1M)„2 + 

K JK. * JK 



1.8.1 Harnack inequalities 



\ (|VRm| 2 + |VRic| 2 ) + |Rm| 2 + |Rzc| 2 



i2\rr... |2 



Vm a | (46b) 



The eigenfunction. We have for 7 = A, u\_ c = \Ju\ + e 2 that 

supuA.e < -D(»7, A)infu A , e (47) 

w w 

The gradient. Now the gradient G\^ e — |Vu A | 2 + e 2 requires that we use the 71 = A + k and we get that 

sup G A . e < D( V , A) (A + D{n, /i)||Rtc||i, w ) inf G A , e (48) 

w w 

The hessian estimate The estimate 

supff A , e < D(r 1 ,X)X 2 (X + D(r 1 ,^\\Ric\\ hW )miH^ (49) 



w 



w 



The estimates for the restriction on the spherical front The restriction of the eigenfunction 

u(r,0) =e A sin(/3i0) 
of the spherical front satifies the following inequalities for j = 1,2: 

f 9 2 ^u\ 2 <C [ 6 2 (X + R 2 )i 
Jf Jf 



)u 2 



where R is a polynomial depending on Rm, VRm, . . . , V 2 Rto, Ric, VRic, VRic. This combined with Michael- 
Simon Sobolev inequality provides Harnack inequalities for the restriction of u, ^Cu\^c 1 u on the spherical front. 

1.8.2 The Bernstein inequalities 

The integration by parts formulas suggest along with the Harnack inequalities the following Berstein type 
inequalities in geodesic pixels: 

Estimates 13. The following estimates hold in a domain inside a geodesic pixel W C P 

|V 2 u A | < sup \W 2 u x \ < C 2 (W)X* +1 (|Vu A | + e) < G 3 (W)A^+ 2 (|« A | + e) (50) 

w 

|V« A | < sup |Vu A | < C 4 (W)At+ 1 (| UA | + e ) (51) 

w 
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